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Abstract 

Background: Prosopis rubriflora and Prosopis ruscifolia are important species in tlie Cliaquenian regions of Brazil. 
Because of the restriction and frequency of their physiognomy, they are excellent models for conservation genetics 
studies. The use of microsatellite markers (Simple Sequence Repeats, SSRs) has become increasingly important in 
recent years and has proven to be a powerful tool for both ecological and molecular studies. 

Findings: In this study, we present the development and characterization of 10 new markers for P. rubriflora and 13 
new markers for P. ruscifolia. The genotyping was performed using 40 P. rubriflora samples and 48 P. ruscifolia 
samples from the Chaquenian remnants in Brazil. The polymorphism information content (PIC) of the P. rubriflora 
markers ranged from 0.073 to 0.791, and no null alleles or deviation from Hardy-Weinberg equilibrium (HW) were 
detected. The PIC values for the P. ruscifolia markers ranged from 0.289 to 0.883, but a departure from HW and null 
alleles were detected for certain loci; however, this departure may have resulted from anthropic activities, such as 
the presence of livestock, which is very common in the remnant areas. 

Conclusions: In this study, we describe novel SSR polymorphic markers that may be helpful in future genetic 
studies of P. rubriflora and P. ruscifolia. 
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Findings 

Background 

The genus Prosopis L. belongs to the Leguminosae 
botanical family, which contains 44 species. Prosopis L. is 
predominantly restricted to the neotropics [1]. Prosopis 
rubriflora [2] and Prosopis ruscifolia [3] are tree species 
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known locally as "espinheiro" and "algarroba," respectively. 
These species are important both economically and eco- 
logically. For example, the fruits and seeds of P, ruscifolia 
are reported to be good sources of nutrition for humans 
and animals [4], and the flowers of P, rubriflora, which are 
present throughout the year, provide important food re- 
sources, such as pollen and nectar, for the local fauna [5]. 
P, rubriflora has a narrow distribution range and is limited 
to Paraguay and Brazil, but P, ruscifolia is also found in 
Argentina and Bolivia [6,7]. 

In Brazil, P, rubriflora and P, ruscifolia are associated 
with Chaquenian areas [8] and are limited to the southern 
portion of the Pantanal [9,10]. Both species are excellent 
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Table 1 Primers developed for Prosopis rubriflora and Prosopis ruscifolia 



Marker GenBank register no. 



Primer sequences (5-3') 



Motifs 



Ta rC) Size^ (bp) Crossed amplification 



Prbl 

Prb2 

PrbS 

Prb4 

PrbS 

Prb6 

Prb7 

PrbS 

Prb9 

PrblO 

Prscl 

Prsc2 

Prsc3 

Prsc4 

PrscS 

Prsc6 

Prsc7 

PrscS 

Prsc9 

Prscl 0 
Prscl 1 



KF923365 

KF923366 

KF923367 

KF92336S 

KF923369 

KF923370 

KF923371 

KF923372 

KF923373 

KF923374 

KC753210 

KC753211 

KC753212 

KC753213 

KC753214 

KC753215 

KC753216 

KC753217 

KC75321S 

KC753219 
KC753220 



F: AACTACCGCAGCAC^CAGA 
R: ACTAOTGGAGATGCCGTGGA 

F: GAAAGCCGCGCTCCTAAG 
R: ATO^GTGTOTGTOTCTCG 
F: TCCAAAGACCGCAAGAAGAT 
R: AGGCCAAAAAGGACTCAAAAT 
F: ATCCGATAAATACACOTCTGG 
R: GGTGTATCGTAAAAGCCTGG 
F: mAAACATOCACGTGAACCTAT 

R: TOACCCCTAAACCCCOT 
F: TCATCTCTCAAAGAAAACGCACTC 
R: CCGCAGAGAAGCCCCTACATA 
F: GGCTOGCATCACCCTCCAT 
R: CTOCCCmCAGTCCAmACCA 
F: CAACACCAAAACGGCGAGATGAT 
R: TOGCCAAACGCCAGCATOG 
F: TTCYTCYCCUCYTCAJCYTCCYCC 
R: ACAACGTOATCCCAAAACCTAAG 

F: ^GGTGGAmGATAGAGCC 
R: GAGTGGGGTCAAGAAAGAACAG 
F: AATGGAGmGmGTGTCTGTGG 
R: ATOCGGATACATCGAGCOTOT 

F: GCGGAATOCAAACGACAA 
R: ACAGCAACACCCTCACTCTCAA 
F: CCACAAGCACACGCACACTCAGAC 
R: CCAGCACTAGAOTCGCCACCAAC 
F: CAAAATCCAACAAATAAACACACC 
R: GGCGGATTOTGGCTCTCT 
F: CGCGTOAGTCTGCOTGCm 
R: CTCATGGTAmCCOTGTCGTCC 
F: CGAGCGGCGAAAAATGATAAA 
R: GCTGOTCCCATAATCCTCTCCT 
F: CAGGGAmAATCTCmGGTGTAG 
R: ACAAGCTGGAAAGAGTCGCA 
F: AGTGACGTGAACACGCTGAGG 
R: TGCTGATGTGTGTGG^GAGAT 
F: TCAGACTCCCGTGAACCAG 
R: CGCACTCGAGCAGCATCT 
F: AACGCAACGGCCGCAACTAT 
R: ACAAAACGCTCGAATACTGGGGG 
F: CCCGGCAACTCAAATCAAOTCATA 



(gt)7 

(gc)4(ac)7 

(ca)7 

(ca)S 

(ac)9 

(tg)lO 

(ac)S 

(gt)13 

(ac)9 

(tca)5 
(ac)9(ct)5 

(ac)9 

(ca)6 



62.7 255-267 



(gt)s 



(gt)s 



61.0 



61.0 



61.0 



564 



56.4 



61.0 



61.0 



62.7 



564 



56.5 



(caa)2(ga)4 63.9 



59.0 



63.9 



(tg)S...(gtgg)2(gt)5 59.0 



(tg)io 
(tg)9 

(ca)7(ct)7 
(ac)ll 



62.7 

59.0 

56.5 
62.7 



140-146 



149-159 



194-230 



149-155 



115-133 



219-225 



144-164 



167-175 



223 



279-297 



64.7 224-252 



64.7 156-160 



21S-232 

220-240 

1S4-210 

122-156 

9S-120 

112-122 

260-2S4 
229-371 



267^ 



126° 



145^ 



203^ 



125^ 



220^ 



154° 



190^ 



300° 



250^ 



122 



230° 



218 



200 



122° 



114° 



244° 
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Table 1 Primers developed for Prosopis rubriflora and Prosopis ruscifolia (Continued) 







R: GGTCTAATOTATOGTGGGCTCTCTGG 










Prscl2 


KC753221 


F: GGGGTGCATGTOGGGATO 
R: mGGCCGGATOAAACAGAGCA 


(gt)io 


59.0 


185-223 


220^ 


PrsclS 


KC753222 


F: OTCACCATCACCGAmCCOT 
R: GCAACGAAGCAGCTGAAGAACAC 


(ctt)5 


62.7 


102 


116 



Ta - Optimal annealing temperature defined after gradient tests of the corresponding markers. ^Range of the fragment sizes from the polymorphic markers and 
the sequenced size of the monomorphic markers; "^Polymorphism observed for the transferred markers based on the 5 samples used. 



indicators of Chaquenian areas in Brazil; P. rubriflora is 
usually associated with arboreal physiognomy, and P. 
ruscifolia is frequently associated with forest physiognomy. 
Both species can be used as models for genetic studies of 
diversity in these areas. 

While estimating genetic diversity, the use of mo- 
lecular markers has been helpful in defining alleles and 
studying genetic flow, population structure, paternity, 
inheritability, genetic maps and conservation genetics 
[11], Simple sequence repeat markers (SSRs), com- 
monly referred to as microsatellite markers, are desir- 
able tools because they are co-dominant in nature, multi- 
allelic and widely distributed in the genome; they are also 
currently cheap, reproducible and relatively easy to analyze 
[12]. This work reports the development, characterization 
and transferability of microsatellite markers for P, rubriflora 
and P, ruscifolia. 

Construction of a microsatellite-enriched library 

DNA was extracted from P, rubriflora and P, ruscifolia 
using the DNeasy® Plant Mini Kit (Qiagen, Hilden, DE) ac- 
cording to the manufacturers instructions. Microsatellite- 
enriched libraries for P, rubriflora and P, ruscifolia were 
constructed as described by Billote et al [13]. The gen- 
omic DNA was digested with Afal after enrichment 



with streptavidin-coated magnetic beads (Streptavidin 
MagneSphere Paramagnetic Particles, Promega, Madison, 
WI); biotinylated (CT)8 and (GT)8 microsatellite 
probes were added for the dinucleotide-enriched li- 
brary. The fragments were amplified by PGR and 
cloned into the pGEM-T vector (Promega, Madison, 
WI). XL 1 -Blue {Escherichia coli) competent cells were 
transformed with the recombinant plasmids and then 
cultivated on agar medium containing ampicillin 
(100 mg/ml), X-galactosidase 2% (100 (ig/ml) and 
IPTG (100 mM). The selected clones were added to a 
Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied 
Biosystems, Foster City, CA) and sequenced using an ABI 
377 sequencer (Applied Biosystems, Foster City, CA). The 
sequences were aligned and edited using SeqMan 
Software (DNAStar, Madison, WI), and the adapters 
and restriction sites were removed using Microsat Software 
(A. M. Risterucci, CIRAD, personal communication). 
To identify microsatellite-enriched regions, we used the 
Simple Sequence Repeat Identification Tool (SSRIT) 
[14] and defined the following numbers of repeats/mo- 
tifs: five/dinucleotides, four/trinucleotides and three/ 
tetra- or pentanucleotides. After these steps, primers 
were designed using the PrimerSelect software (DNAS- 
tar, Madison, WI). 



Table 2 Markers developed for Prosopis rubriflora 



Marker 




Na 






Ho 




He 


PIC 


Null alleles 


HW (P-value) 


T 


FSM 


FSV 


FSM 


FSV 


FSM 


FSV 


FSM 


FSV 


FSM 


FSV 


Prbl 


7 


6 


7 


0.526' 


0.474 


0.522' 


0.737 


0.602 


0.044 


0.145 


0.348 


0.045 


Prb2 


3 


1 


3 


0.000 


0.150' 


0.000 


0.145' 


0.073 


0.001 


0.000 




1.000 


Prb3 


3 


3 


2 


0.316 


0.150 


0.428 


0.296 


0.313 


0.103 


0.134 


0.101 


0.069 


Prb4 


12 


11 


8 


0.650 


0.850' 


0.799 


0.803' 


0.766 


0.092 


0.000 


0.457 


0.291 


Prb5 


4 


4 


4 


0.500 


0.600 


0.583 


0.683 


0.576 


0.039 


0.007 


0.689 


0.254 


Prb6 


6 


4 


6 


0.500' 


0.400 


0.458' 


0.432 


0.413 


0.000 


0.000 


1.000 


0.335 


Prb7 


4 


4 


4 


0.684' 


0.650' 


0.605' 


0.499' 


0.473 


0.000 


0.000 


0.355 


0.121 


Prb8 


5 


5 


3 


0.500 


0.250 


0.524 


0.304 


0.370 


0.000 


0.050 


0.821 


0.468 


Prb9 


10 


10 


9 


0.684 


0.650 


0.835 


0.819 


0.791 


0.053 


0.062 


0.177 


0.016 



FSM - Fazenda Sao Manoel, FSV - Fazenda Santa Vergmia, Ng - Number of alleles, Hq - Observed heterozygosity, Hg - Expected heterozygosity, PIC - Polymorphism 
information content, P-values of Hardy-Weinberg (HW) equilibrium (P-value > 0.0055 after Bonferroni correction), null alleles (null frequency < 0.20). ^Populations 
where the values of Hq were higher than those of He. 
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Table 3 Markers developed for Prosopis ruscifolia 



Marker 




N 




H 




H 




PIC 


Null alleles 


HW (P-value) 


T 


FRC 


ECD 


FRC 


ECD 


FRC 


ECD 


FRC 


ECD 


FRC 


ECD 


Prscl 


7 


4 


7 


0.600' 


0.783 


0.566' 


0.828 


0.701 


0.000 


0.020 


0.916 


0.009 


Prsc2 


8 


7 


5 


0.480 


0.565 


0.553 


0.761 


0.658 


0.032 


0.114 


0.034 


0.00' 


PrscS 


3 


3 


3 


0.360' 


0.348 


0.344' 


0.456 


0.348 


0.000 


0.048 


0.674 


0.073 


Prsc4 


5 


4 


3 


0.040 


0.043 


0.321 


0.275 


0.289 


0.247^ 


0.221^ 


0.000' 


0.000' 


Drcr-Q 

r 1 SCj 


0 


4 


r 
J 


U.44U 


U.z/o 


U.j/o 


U.Dj4 


n A Q/i 
U.4o4 


0.000 


0.218^ 


1.000 


0,000' 


Prsc6 


8 


4 


8 


0.320 


0.696 


0.653 


0.801 


0.703 


0.196 


0.068 


0.001' 


0.087 


Prsc7 


9 


6 


7 


0.440 


0.565 


0.727 


0.779 


0.755 


0.167 


0.081 


0.001' 


0.000' 


PrscS 


7 


4 


7 


0.320 


0.174 


0.577 


0.789 


0.656 


0.165 


0.338^ 


0.008 


0.000' 


Prsc9 


5 


3 


3 


0.240 


0.130 


0.280 


0.559 


0.430 


0.040 


0.267^ 


0.484 


0.000' 


Prscl 0 


7 


5 


5 


0.286 


0.227 


0.633 


0.758 


0.670 


0.212^ 


0.294^ 


0.000' 


0.000' 


Prscl 1 


17 


8 


12 


0.560 


0.652 


0.845 


0.884 


0.883 


0.157 


0.127 


0.003 


0.000' 


Prscl 2 


5 


4 


5 


0.600 


0.591 


0.569 


0.707 


0.589 


0.000 


0.075 


0.592 


0.013 



FRC - Fazenda Retire Concei^ao, ECD - Esta^ao do Carandazal, Na - Number of alleles, Hq - Observed heterozygosity. He - Expected heterozygosity, PIC - Polymorph- 
ism information content, P-values of Hardy-Weinberg (HW) equilibrium (P-value > 0.0041 after Bonferroni correction). ^Populations where the values of Hq were 
higher than those of Hg; "^Possible null alleles (null frequency < 0.20); "^Departure from HW equilibrium was observed. 



Fragment amplification 

The fragments were amplified using polymerase chain 
reactions containing 8 ng of template DNA, 2 mM MgCl2, 
50 mM KCl, 20 mM Tris-HCl (pH 8.4), 0.2 mM dNTPs, 
0.19 mg/ml BSA (bovine serum albumin), 0.15 mM of 
each primer and 1 U of Taq DNA polymerase; the reac- 
tions were then brought to a final volume of 20 (il with 
ultrapure water. To define the temperatures for the PGR 
reactions, we adopted the guidelines described by 
Mottura et al, [15]; for the annealing temperatures, we 
used a gradient program with temperatures ranging from 
65°C to 55°C. The samples were collected in the Chaco 
remnants of Corumba and Porto Murtinho, Mato Grosso 
do Sul, Brazil. Twenty P, rubriflora samples were collected 
in each of two Chaco remnant locations: Fazenda Sao Man- 
oel (FSM) (21°47'44.5"S; 57°39'34.6"W) and Fazenda 
Santa Verginia (FSV) (22°06'40.5"S; 57°49'57.6"W). 
Twenty-three P, ruscifolia samples were collected in 
Esta^ao do Carandazal (ECD) (19°48'33.2"S; 57°10'11.0" 
W), and 25 samples were collected in Fazenda Retiro 
Concei^ao (FRC) (21°42'23.7"S; 57°45'58.2"W). The 
cross-amplification of the markers was evaluated in 5 P, 
rubriflora samples obtained from FRC (21°41 '00.7"S; 
57°46'43.8"W) and 5 P, ruscifolia samples from 
Chacara Jacare (21°41 '20.1 "S; 57°52'15.5"W) using 
the same conditions as for the native species. The amp- 
lified samples were genotyped by vertical electrophor- 
esis using denaturating polyacrylamide gels (6%), and 
DNA bands were visualized using silver nitrate [16]; the 
sizes of the resulting fragments were estimated by compari- 
son to a 10-bp DNA ladder (Invitrogen, Carlsbad, CA). Stat- 
istical analyses were performed using Microsatellite Toolkit 
v.3.1.1 [17] to calculate the expected heterozygosity (Hg), 



observed heterozygosity (Ho) and polymorphism infor- 
mation content (PIC). The Genepop software v. 1.2 [18] 
was used to estimate adherence to Hardy-Weinberg 
(HW) equilibrium and possible linkage disequilibrium 
(LD), and the frequency of null alleles was estimated using 
FreeNA [19]. 

Results and discussion 

We designed 32 primer pairs: 13 for P, rubriflora and 19 
for P, ruscifolia. However, only 10 of the P, rubriflora 
primer pairs and 13 of the P, ruscifolia primer pairs 
amplified properly. The nine remaining pairs of primers 
were discarded because amplification errors were observed 
in the preliminary tests. Polymorphisms were detected in 9 
of the native P, rubriflora markers and 12 of the native P, 
ruscifolia markers; only one marker from each species had 
a monomorphic pattern based on the populations analyzed. 
Eight markers from P. rubriflora successfully cross- 
amplified and were polymorphic for the tested samples, 
and 2 markers failed during cross-amplification. Eleven P, 
ruscifolia markers were successfully cross-amplified; 7 
were polymorphic, and 2 failed this analysis (Table 1). 

The number of P. rubriflora alleles in the sampled rem- 
nants ranged from 3 to 12; the polymorphism information 
content (PIC) values of these markers ranged from 
0.073 to 0.791, the observed heterozygosity (Ho) ranged 
from 0.000 to 0.850, and the expected heterozygosity (Hg) 
ranged from 0.000 to 0.835. No evidence of null alleles 
was observed, and no departure from Hardy-Weinberg 
equilibrium was observed (Table 2). No significant link- 
age disequilibrium (LD) was observed for any of the 
markers of this species after Bonferroni correction 
(P-value for 5% = 0.001389). The number of P. ruscifolia 
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alleles in both of the remnants ranged from 3 to 17, the 
PIC values ranged from 0.289 to 0.883, the Ho values 
ranged from 0.040 to 0.783, and the Hg values ranged 
from 0.275 to 0.884. Possible null alleles were observed 
for the markers PrscS, Prsc8 and Prsc9 from one remnant 
(ECD), and the markers Prsc4 and PrsclO had possible 
null alleles in both remnants. A departure from HW was 
observed for Prsc2, PrscS, Prsc6, Prsc8, Prsc9 and Prscll 
in one of the remnants (the majority were observed in 
ECD) and for Prsc4, Prsc7 and PrsclO in both remnants 
(Table 3). Significant LD was observed for the loci 
PrscS and Prsc6 after Bonferroni correction (P-value for 
1% = 0.00016). 

Higher values of Ho were observed for the Prbl, Prb2, 
Prb4, Prb6, Prb7, Prscl, Prsc3 and PrscS markers in this 
study; these higher values may indicate that an insuffi- 
cient number of samples was collected or may be related 
to the reproductive patterns of these populations. The 
ECD populations are highly disturbed, and the PRC 
population is currently recovering from a relatively re- 
cent suppression (within the last IS years); these fac- 
tors may underlie the observed departure from HW 
and the presence of null alleles. A study with new and 
conserved populations may produce better results for 
these markers. 

These markers are the first microsatellite markers devel- 
oped for Prosopis rubriflora and Prosopis ruscifolia, and 
together with the set of P, ruscifolia markers amplified by 
Bessega et al [20], they are expected to be useful tools for 
studies of the conservation genetics, reproductive biology, 
phylogeography and taxonomy of these species. 

Availability of supporting data 

The original sequences of the developed markers were 
submitted to the GenBank database (http://ncbi.nlm.nih. 
gov), and the registered codes are available in Table 1. 

The testimony samples were deposited at Herbarium 
Universidade Estadual de Campinas (UEC - Campinas, SP, 
BR) and registered according to the following: P, rubri- 
flora - l^^l (Fazenda Sao Manoel - Porto Murtinho, 
MS), 1S471S (Fazenda Santa Verginia - Porto Murtinho, 
MS); P, ruscifolia - 74469 (Fazenda Retiro Concei^ao - 
Porto Murtinho, MS), 37266 (Esta^ao do Carandazal - 
Corumba, MS). 
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